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The Binary-Encounter-Bethe (BEB) model 
for electron-impact total ionization cross 
sections has been applied to CHJ, CH3, 
CHJ, C2H2, C2HJ, C2HJ, and H3O*. The 
cross sections for the hydrocarbon ions are 
needed for modeling cool plasmas in fu- 
sion devices. No experimental data are 
available for direct comparison. Molecu- 
lar constants to generate total ionization 
cross sections at arbitrary incident elec- 
tron energies using the BEB formula are 
presented. A recent experimental result 
on the ionization of H3O* is found to be al- 
most 1/20 of the present theory at the 
cross section peak. 
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1. Introduction 

Ionization cross sections for atomic and molecular 
ions are among the critical data needed in modeling 
plasmas in fusion devices. Hydrocarbon molecules and 
their ion fragments are formed inside a tokamak in edge 
plasmas and near a divertor. The Binary-Encounter- 
Bethe (BEB) model [1] has successfully generated reli- 
able total ionization cross sections of small as well as 
large molecules [2-6]. The BEB model combines a mod- 
ified form of the Mott cross section with the asymptotic 
form of the Bethe theory (i.e., high incident energy T) 
for electron-impact ionization of a neutral atom or 
molecule. The original BEB model was slightly modi- 
fied for applications to atomic and molecular ions [7]. 

In this article we apply the modified BEB formula for 
ions to hydrocarbon ions of interest to magnetic fusion: 
CHJ, CH^, CH:, C2HJ, C2H;, C2HJ, and HjO^ We out- 
line the theory in Sec. 2, and our theoretical results are 
presented in Sec. 3. A recent experiment on the forma- 



tion of H3O** by electron impact [8] is compared to the 
present theory in Sec. 3. 



2, Outline of Theory 

The BEB formula for ionizing an electron from a 
molecular orbital of a neutral molecule by electron im- 
pact is [1]: 



Cbeb ■ 



t+u+\ 



ln_f 



+ 1 



Inf 



t+ 1 



(1) 



where t=T/B, u=U/B, S = 4TTalN RVB'-, a,, is the 
Bohr radius (= 0.5292 A), R is the Rydberg energy 
(= 13.6057 eV), Tis the incident electron energy, and A^, 
B , and U are the electron occupation number, the bind- 
ing energy, and the average kinetic energy of the orbital, 
respectively. 
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In Eq. (1), the terms in the square brackets are based 
on the Mott theory and the Bethe theory. However, the 
denominator t + u + I is based on a plausible, but less 
rigorous argument, i.e., the effective kinetic energy of 
the incident electron seen by the bound target electron 
should be the incident electron energy T plus the poten- 
tial energy U + B oi the target electron [9]. Hence the T 
in the denominator of the original Mott and Bethe theo- 
ries was replaced hyT+U + B,ort + u+lin Eq. (1), 
where B is used as the energy unit. 

The net effect of using t+u+\ instead of t in the 
denominator of Eq. (1) is to reduce substantially the 
cross section near the ionization threshold. This modifi- 
cation was found not only to be effective but also abso- 
lutely necessary to have the theory agree with reliable 
experimental ionization cross sections near the threshold 
for many neutral atoms and molecules. 

In a previous article [7] for singly charged molecular 
ions, we have shown that the denominator t + u+\ is 
replaced by t + {u+ l)/2 to generate ionization cross 
sections in good agreement with available experimental 
data. The modified BEB equation for singly charged 



t+{u+ l)/2 



Inf 



+ 1 



Inf 

7TT 



(2) 



Equation (2) is as simple as the BEB formula for neutral 
targets, Eq. (1), and does not require any more input 
data than the original BEB formula. 



3. Theoretical Results 

We present the BEB cross sections from Eq. (2) for 
cm, cm, CHJ, C2HJ, CjHt CzHJ, and H3O* in Figs. 
1-4. The molecular constants B, U, and A^ for the 
molecules are listed in Table 1. For all molecular ions 
except CHJ, molecular geometries were computed using 
a hybrid density functional (B3LYP) [10,11] with 6- 
31G(d) basis sets. For CHJ, B3LYP/6-31G(d) gave an 
incorrect molecular symmetry (C2 point group instead 
of C2v), so the geometry was computed using frozen- 
core, second-order perturbation (MP2) theory with 6- 
31G(d) basis sets. The B3LYP or MP2 geometries were 
used for all subsequent calculations of B and U. Kinetic 
energies U for all orbitals, and binding energies B for the 
inner orbitals, were calculated at the Hartree-Fock (HF) 
level using 6-311G(d,p) basis sets. More accurate, cor- 
related values of B were obtained for the outer-valence 
orbitals by using frozen-core Green's function (OVGF) 
methods [12,13] and 6-311-i-G(d,p) basis sets. For the 
important threshold ionization, B values were obtained 
by using frozen-core coupled cluster theory [CCSD(T)], 
with the single and double excitation operators included 
iteratively [14] and the contribution from connected 
triples estimated perturbatively [15]. Dunning's correla- 
tion-consistent valence-triple-zeta (cc-pVTZ) basis sets 
[16] were used for the CCSD(T) calculations. The HF 
calculations were performed using the GAMESS [17] 
program package; all other calculations employed Gaus- 
sian 98' [18]. 




Fig. 1. Electron-impact ionization cross sections for CHJ and C2HJ. 



Certain commercial equipment, instruments, or materials are identified 
in tills paper to foster understanding. Sucii identification does not imply 
recommendation or endorsement by the National Institute of Standards 
and Technology, nor does it imply that the materials or equipment iden- 
tified are necessarily the best available for the purpose. 
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Fig. 2. Electron-impact ionization cross sections for CH* and CzHj, 
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Fig. 3. Electron-impact ionization cross section for CH4 and C2HJ. 



•< 




Fig. 4. Electron-impact ionization cross section for H3O*. 
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Table 1. Molecular point group, molecular orbitals (MO), electron 
binding energy B, kinetic energy U , and electron occupation number 
N for CHJ, CH3, CHJ, C2HJ, CiHj, C2HJ, and H3O* 



Molecule (point group) 



MO 



B(eV) 



CHt (C2,,) 



CHJ (D3,,) 



CH: (C2..) 



C2HJ (D.a) 



C2H: (D2) 



C2HJ (D3,) 



H3O* (C3„) 



C/(eV) 



TV 



la, 


319.07 


436.58 


2 


2fl, 


33.16 


38.98 


2 


1^2 


27.04 


29.90 


2 


3a, 


22.17 


37.74 


1 


\a[ 


317.93 


436.46 


2 


la\ 


33.55 


37.57 


2 


le' 


25.59 


29.56 


4 


la, 


315.92 


436.20 


2 


2a, 


33.67 


33.40 


2 


3a, 


25.87 


27.79 


2 


\b2 


24.48 


28.73 


2 


lb. 


22.08 


29.42 


1 


la. 


316.56 


435.87 


2 


la-„ 


316.48 


436.63 


2 


2a-. 


34.59 


49.00 


2 


2o-„ 


27.96 


35.66 


2 


3a-, 


26.06 


34.23 


2 


l77„ 


21.00 


31.42 


3 


la 


315.15 


436.08 


2 


lb. 


315.12 


436.44 


2 


la 


33.52 


40.73 


2 


lb. 


28.22 


36.01 


2 


lb2 


24.41 


27.06 


2 


3a 


23.30 


35.07 


2 


1^3 


21.81 


29.99 


2 


2^3 


19.21 


29.54 


1 


Iflig 


314.28 


436.27 


2 


Iflju 


314.28 


436.30 


2 


2a„ 


31.76 


33.19 


2 


2a2u 


29.16 


38.10 


2 


le„ 


23.07 


26.29 


4 


1., 


21.38 


30.07 


4 


3aig 


19.04 


30.07 


1 


la. 


571.51 


794.48 


2 


2a, 


45.59 


72.42 


2 


le 


29.90 


52.11 


4 


3a, 


24.7" 


65.41 


2 



" Experimental value from Ref. [8]. 

In general, when an electron collides with a molecu- 
lar ion we get 



e" + AB* ^ A + B* + e , 
or A* + B + e". 

e" + AB* -^ AB** + 2e-, 

or A* + B* + 2e", 



(3) 
(4) 
(5) 
(6) 



or A** + B + 2e" 



or A + B** + 2e 



(7) 



(8) 



Processes (3) and (4) are dissociation without ioniza- 
tion, while processes (5) through (8) are the ionizing 
events described by the BEB model. The model calcu- 
lates the sum of all processes (5) through (8) that lead 
to the ejection of a bound electron. Moreover, the model 
also assumes — erroneously — that all energy transfers 
from the incident electron to the target molecule that 
exceed the ionization energy of a given molecular or- 
bital result in ionization. This is an assumption common 
to all binary-encounter type theories. Although such an 
assumption may be valid for atoms, molecules may 
dissociate without ionizing even if energy transfers ex- 
ceed the orbital binding energies. If processes (3) and 
(4) are significant for energy transfers above orbital 
binding energies, then the BEB model will overestimate 
ionization cross sections. More discussions on this point 
can be found in Ref. [5]. 

Experimentally, the production of doubly charged 
ions can be detected directly when the doubly charged 
ions have reasonably long lifetimes. In reality, most 
doubly charged molecular ions quickly dissociate into 
two singly charged fragments, making it almost impos- 
sible to distinguish processes (3) and (4) from the disso- 
ciation of doubly charged ions by Coulomb repulsion, 
i.e., process (6). For this reason, it is difficult to distin- 
guish processes (3) and (4) from process (6) simply by 
detecting singly charged ions unless coincidence mea- 
surements of all products are performed. The usual ex- 
perimental procedure is to measure the cross section for 
producing any ion, i.e., (3) through (8). Then, processes 
(3) and (4) are measured separately, and subtracted 
from the total ion production cross section. This sub- 
traction introduces large uncertainties in the resulting 
experimental ionization cross sections. 

The direct measurement of only the doubly charged 
ions — processes (5), (7), and (8) — tends to produce 
small cross sections compared to the total ionization 
cross section because of the high probability for the 
rapid break-up of the doubly charged ions as shown in 
Ref. [7] for CO*. As another example, Bahati et al. [8] 
recently reported experimental cross sections for the 
process of 



e" + H3O* -^ H3O** + le." 



(9) 



Their measurement corresponds to process (5) only, and 
their peak cross section is 0.049 A^ at T ~ 125 eV. The 
position of the peak is in agreement with the BEB cross 
section in Fig. 4, but the magnitude is almost 1/20 of the 
BEB cross section. This discrepancy and the similar 
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discrepancy in CO* (experiment is lower by a factor of 
1/12 at the peak) measured by the same group [19] is a 
strong indication that either dissociative ionization is the 
dominant process, or the break-up of the doubly charged 
ions is faster than the experimental capability to detect. 
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